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ABSTRACT 

CH 4 is proposed to be the starting point of a rich organic chemistry. Solid 
CH4 abundances have previously been determined mostly toward high mass star 
forming regions. Spitzer/IRS now provides a unique opportunity to probe solid 
CH 4 toward low mass star forming regions as well. Infrared spectra from the 
Spitzer Space Telescope are presented to determine the solid CH 4 abundance to- 
ward a large sample of low mass young stellar objects. 25 out of 52 ice sources in 
the did (cores to disks) legacy have an absorption feature at 7.7 /im, attributed 
to the bending mode of solid CH 4 . The solid CH 4 / H 2 abundances are 2-8%, 
except for three sources with abundances as high as 11-13%. These latter sources 
have relatively large uncertainties due to small total ice column densities. Toward 
sources with H 2 column densities above 2xl0 18 cm -2 , the CH 4 abundances (20 
out of 25) are nearly constant at 4.7±1.6%. Correlation plots with solid H2O, 
CH3OH, CO2 and CO column densities and abundances relative to H2O reveal 
a closer relationship of solid CH 4 with CO2 and H 2 than with solid CO and 
CH3OH. The inferred solid CH 4 abundances are consistent with models where 
CH 4 is formed through sequential hydrogenation of C on grain surfaces. Finally 
the equal or higher abundances toward low mass young stellar objects compared 
with high mass objects and the correlation studies support this formation path- 
way as well, but not the two competing theories: formation from CH 3 OH and 
formation in gas phase with subsequent freeze-out. 
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Introduction 



The presence and origin of complex organic molecules in protostellar regions and their 
possible incorporation in protoplanetary disks is an active topic of research. CH 4 is pro- 
posed to be a start ing point of a rich chemistry, especially when UV photons are present 
( IDartois et al.ll2005l ). I n particular CH4 is bel ieved to play a key role in the formation process 
of prebiotic molecules (IMarkwick et al.ll2000l ). 



CH4 is less well studied in interstellar and circumstellar media compared to other small 
organic molecules because CH 4 has no permanent dipole moment and therefore cannot be 
observed by pure rotational transitions at radio wavelengt hs. Solid CH 4 wa s first detected 
through its bending mode at 7.67 fim from the ground by lLacy et al.l (1199 if ) , and with the 



Infrar ed Space Observatory Short Wavelength Spectrometer (ISO-SWS) by iBoogert et al. 



( 119961 ) toward a few high-mass sources. Tentative claims have been made toward some other 
objects incl uding low-mass protost a rs, but are inconclu s ive because of the low S/N ratio of 
these data jCernicharo et al.l I2OO0I : loiirtler et al.l I2OO2I : klexander et alil2003l ). Solid CH 4 
has also been detected from the ground through its stretching mode at 3.3 fim, but only 
toward the brightest high m ass sources due to pr oblems in removing the many atmospheric 
lines in this spectral region (IBoogert et al.ll2004al ). 



Models predict CH 4 to form rapidly on cool grains through successi ve hydrogenation o 



atomic C; similarly H 2 Q is formed through hydrogenation of atomic O ( 



Allen & Robinson! 


1977 




Tielens & Haeen 


Aikawa et al. 


2005 


). Observations of CH 



19821 : iBrown et al.l Il988t lHasegawa et al 



van de Hulst 



1946 



1992 



of grain surface chemistry. Compared to H 2 the observed gas- and solid-state CH 4 abun- 
dances are low; reported CH 4 abunda nces are typically a few percent with respect to H 2 
( lLacy et al.lll99ll ; IBoogert et al.lll998l ). This points to relatively low atomic C abundances at 
the time of CH 4 formation, with most C already locked up in CO as H readily reacts with C 
on surfaces (IHiraoka et al.lll998l ). This is in agreeme nt with the high CH3OH abundances in 



sever al lines of sight, formed by hydrogenation of CO (IDartois et al.lll999l : lPontoppidan et al. 



20031 ). and large C0 2 abundances, formed through oxidation of CO or hydrogenated CO. 
That these molecules are all formed through a similar process is corroborated by the profiles 
of solid C0 2 absorption bands, which usually show an intimate mix t ure of CO2, CH3O H and 
H 2 in interstellar ices (IGerakines et al.lll999l : IBoogert et al.ll2000t iKnez et al.ll2005l ). 
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If CH4 is formed efficiently through grain surface reactions as well, CH4 should be 
similarly mixed with H 2 0. Observations of solid CH 4 toward a few high mass young stellar 
objects (YSOs) show that the CH 4 absorption band profiles are broad and agree better 
with CH 4 in a h ydrogen bonding ic e, H 2 or CH 3 OH, than with a pure CH 4 ice or CH 4 
mixed with CO (IBoogert et al.lll997l ). This profile analysis does not , however, exclude CH 4 
formation from photoprocessing of CH3OH (lAllamandola et al.lll988l ; iGerakines et al.lll996l ). 
In addition, because of the small sample in previous studies, it is unclear if these broad profiles 
are a universal feature. Hence it cannot be excluded that CH 4 in some environments may 
form in the gas phase and subsequently freeze out. 

Because formation pathway efficiency depends on environment, another method for 
testing formation routes is through exploring the distribution of CH 4 toward a large sample 
of objects of different ages, luminosities and ice column densities. In addition, correlations, 
or lack thereof, with other ice constituents may provide important clues to how the molecule 
is formed. If CH 4 is formed through hydrogenation on grain surfaces in quiescent clouds, 
the CH 4 abundance with respect to H 2 should be fairly source independent since this 
mechanism mainly depends on the initial conditions of the cloud before the star forms, 
which seem to vary little between different star forming regions. Because of the generality 
of this mechanism, the CH 4 and H 2 0, and possibly C0 2 , column densities should correlate 
over a large range of different environments. This is also the prediction of several models 
where t he solid CH 4 /H 2 and CH 4 /C0 2 ratios in dark clouds vary little both as a function 



of tim e (IHasegawa et al.lll992l ) and distance into a cloud that has collapsed (lAikawa et al. 
20051 ). In contrast, CO, which is formed in the gas phase and subsequently frozen out, is 



predicted to only correlate with CH 4 during certain time intervals. If CH 4 instead forms in 
the gas phase, solid CH 4 should be better correlated with solid CO than with solid H 2 and 
C0 2 , since pure CH 4 freezes out and de sorbs at similar temperatures to CO (Fraser et al. 
submitted to ASzA, ICollings et al.l I2004J ) . Finally, if CH 4 forms by UV photoprocessing of 
CH3OH more CH 4 would also be expected to form toward sources with stronger UV- fields 
i.e. higher mass objects. 

The objective of this study is to determine the CH 4 abundances and distribution pattern 
toward a sample of low mass young stellar objects, varying in evolutionary stage and total 
ice column density. The distribution pattern and correlations with other ice constituents 
within the sample as well as comparison with high mass young stellar objects will be used to 
constrain the CH 4 formation mechanism. This study is based on spectra acquired with the 
Spitzer Infrared Spectrometer (IRS) as part of our legacy program 'From molecular cores to 
protoplanetary disks' (c2d), whic h provides a large s ample (41 sources) of infrared spectra of 
low mass star formation regions (jEvans et al.ll2003l ). In addition If sources are added from 
the GTO program 2 for which ground based 3-5 /mi already exist (IPontoppidan et al.ll2003l ). 
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Overviews of the H2O, CO2, CH3OH and other ice species in these data are found in Boogert 
et al. (ApJ submitted, Paper I) and Pontoppidan et al. (ApJ submitted, Paper II). The 



detectio n of solid CH 4 toward one of the sources, HH46 IRS, was published by lBoogert et al. 



( I2004bl ). We have detected an absorption feature, which is attributed to solid CH 4 , toward 
25 out of 52 low mass ice sources found in this did sample. 



2. Source sample selection, observations and data reduction 

The source sample consists of a combination of known low mass protostars and new 
protostars discovered by their Spitzer IRAC and MIPS broad-band SEDs. Spitzer/IRS 
spectra were obtained as part of the did Legacy program (PIDs 172 and 179) as well as a 
dedicated open time program (PID 20604) and a few archival spectra observed as a part of 
the GTO programs of Houck et al. Among the targets observed with the IRS short-long 
(SL) module, 41 ice sources were identified from their spectra in the did sample (Paper I). 
The GTO sources are all associated with the Ophiuchus clo ud and were selected based on 



previous ice observations at 3-5/im (IPontoppidan et al.ll2003l ). The source sample of 25 low 



mass protostars presented here was selected from the 41 did and 11 GTO ice sources based 
solely on the existence of an absorption feature at 7.7 /zm, identified with solid CH 4 (Table 
[TJ. Due to the high sensitivity of the Spitzer/IRS, a large range of star formation stages 
are represented in the sample from very young YSOs at the interface of the Class 0/1 stages 
(e.g. Bl-c) to objects like RNO 91, identified as an inclined disk with a remn ant envelope. 



It als o includes VEry Low Luminosity Objects (VELLOs) such as L1014 IRS (jYoung et al. 



20041 ) . More information about the evolutionary stages of the objects are reported in Paper 



I. 

The Spitzer/IRS spectra were reduced using the 2-dimensional Basic Calibrated Data 
(BCD) spectral images, produced by the Spitzer Science Center (SSC) pipeline version 
S13.2.0, as a starting point. The low resolution modules, which are the relevant ones for the 
solid CH 4 feature, were reduced in a similar way to ground-based spectra. First, the spectral 
and spatial dimensions were orthogonalized, and then the 2-dimensional images of the two 
nodding positions were subtracted in order to remove the extended emission. A fixed width 
extraction was performed and then the 1-dimensional spectra were averaged. Subsequently 
the spectra were divided by spectra of the standard star HR 2194, reduced in the same way, 
to correct for wavelength dependent slit losses. Finally the spectra were multiplied along the 
flux scale in order to match Spitzer/IRAC photometry. One source, EC 92, was carefully 
extracted manually due to multiple other sources present in the slit. More details about the 
reduction and the complete mid-infrared spectra of the did sources is presented in Paper I. 
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The CH 4 7.7 /im absorption feature falls within the range of the IRS modules SL1 and SL2. 
At these wavelengths the resolving power R = A/AA of IRS is 65 or 0.12 /im for SL1 and 
125 or 0.06 /im for SL2 (IRS Data Handbook). The plotted spectra are made up using data 
from the main orders of SL1 and SL2, while the SL1 bonus order was only used to confirm 
detections. SL2 stops around 7.6 /im so the peak of the CH4 absorption feature is always in 
SL1. The part of the spectra from the SL2 module is hence mainly used to determine the 
continuum and in some cases the shape of a low wavelength wing. When comparing obser- 
vations to laboratory spectra, the laboratory spectra are always convolved to the resolution 
of SL1 (0.12 /xm, sampled with two pixels per resolution element). 

ISO-SWS spectra of four high mass YSOs are used for comparison between low mass 
and high mass protostars in this paper. The latest Standard Processed Data (SPD) pipeline 
version 10.1 products were taken from the ISO archive and the detector scans were cleaned 
from cosmic ray hits and averaged. The final spectra do not show significant differences with 
respect to the data published in Keane et al. (2001). 
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Table 1. The source sample of 25 low mass stars observed with Spitzer-IRS. In addition 4 
high mass stars with previously published ISO spectra and CH 4 detections are included for 

comparison. 



Source 


Alias 


RA J2000 


Doc J2000 


Cloud 


Type 


Obs. ID 


IHAb 03235+3004 




no o/" 1 ot rr 

03:2b:37.5 


i on 1 et ot n 

+30:15:27.9 


Perseus 


low 


no o er cr on 

9835520 


IKAb 03245+3002 




03:2 ( :39.0 


+30:12:59.3 


Perseus 


low 


6368U00 


T i jrc CT\ ttv T-i 

1j1400 blVilVli 




03:2 ( :43.3 


i on. 1 o . oo q 
+30:12:28.8 


Perseus 


low 


1591 (056 


IRAS 03254+3050 




03:28:34.2 


+31:00:51.2 


Perseus 


low 


11827200 


Bl-c 




03:33:17.9 


+31:09:31:0 


Perseus 


low 


13460480 


Bl-b 




03:33:20.3 


+31:07:21.4 


Perseus 


low 


1596544 


L1489 IRS 


IRAS 04016+2610 


04:04:43.1 


+26:18:56.4 


Taurus 


low 


3528960 


IRAS 08242-5050 


HH46 IRS 


08:25:43.8 


-51:00:35.6 


HH46 


low 


5638912 


IRAS 12553-7651 




12:59:06.6 


-77:07:40.0 


Cha 


low 


9830912 


IRAS 15398-3359 




15:43:02.3 


-34:09:06.7 


B228 


low 


5828864 


GSS 30 IRS1 




16:26:21.4 


-24:23:04.1 


Ophiuclms 


low 


12699392 


IRS 42 




16:27:21.5 


-24:41:43.1 


Ophiuchus 


low 


12699648 


IRS 43 




16:27:27.0 


-24:40:52.0 


Ophiuclms 


low 


12699648 


IRS 44 




16:27:28.1 


-24:39:35.0 


Ophiuclms 


low 


12699648 


IRS 63 




16:31:35.7 


-24:01:29.5 


Ophiuchus 


low 


12676608 


VSSG 17 


IRS 47 


16:27:30.2 


-24:27:43.4 


Ophiuchus 


low 


12698624 


RNO 91 


IRAS 16316-1540 


16:34:29.3 


-15:47:01.4 


L43 


low 


5650432 


B59 YS05 




17:11:22.2 


-27:26:02.3 


B59 


low 


14894336 


2MASSJ17112317-2724315 




17:11:23.1 


-27:24:32.6 


B59 


low 


14894592 


SVS 4-5 


EC 88 


18:29:57.6 


+01:13:00.6 


Serpens 


low 


9407232 


EC 92 


SVS 4-10 


18:29:57.9 


+0.1:12:51.6 


Serpens 


low 


9407232 


CrA IRS 5 




19:01:48.0 


-36:57:21.6 


Corona Australis 


low 


9835264 


CrA IRS 7 B 




19:01:56.4 


-36:57:28.0 


Corona Australis 


low 


9835008 


CrA IRAS32 




19:02:58.7 


-37:07:34.5 


Corona Australis 


low 


9832192 


L1014 IRS 




21:24:07.5 


+49:59:09.0 


L1014 


low 


12116736 


W33A 




18:14:39.4 


-17:52:01.3 




high 


32900920 


NGC7538 IRS9 




23:14:01.6 


+61:27:20.2 




high 


09801532 


GL989 




06:41:10.1 


+0.9:29:35.8 




high 


71602619 


GL7009S 




18:34:20.9 


-05:59:42.2 




high 


15201140 
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Table 2. Gaussian parameters of the observed absorption features and the CH 4 column 
densities and abundances relative to solid H 2 and upper limits of S0 2 .. 



Source 


A 


FWHM 


T pcak 


Jt 7 . 7 


/ T CH 4 


N(CH 4 ) N(CH 4 ) N(S0 2 )ma X 




/an 






cm -1 


cm -1 


10 17 cm- 2 /N(H 2 O)xl00 



IRAS 03235+3004 


7. 


.69+0 


,03 





16+0, 


05 


0.16+0.07 


4.6+1. 





2.9+0 


,8 


6.2+1.8 


4.3+1.4 


0.35+0.19 


IRAS 03245+3002 


7. 


.59+0 


03 





.20+0, 


08 


0.17+0.04 


4.3+0. 


1 


3.2+0 


,4 


6.8+0.8 


1.7+0.3 


0.08+0.01 


L1455 SMM1 


7. 


.66+0. 


01 


0. 


15+0. 


02 


0.24+0.03 


6.2+0. 


3 


4.9+0. 


,1 


10.+0.1 


5.8+0.9 


0.21+0.03 


IRAS 03254+3050 


7. 


.68+0 


,01 





20+0. 


03 


0.040+0.008 


1.4+0. 


2 


0.7+0 


,1 


1.6+0.3 


4.0+0.9 


0.45+0.14 


Bl-c 


7. 


.68+0 


.01 


0, 


13+0. 


01 


0.40+0.05 


8.7+1. 


3 


7.6+1 


,2 


16+3 


5.4+1.4 


0.11+0.03 


Bl-b 


7. 


.70+0. 


.01 


0, 


10+0. 


01 


0.16+0.01 


2.7+0. 


1 


2.8+0. 


,1 


5.9+0.2 


3.3+0.6 


0.00+0.01 


L1489 IRS 


7. 


.68+0 


01 


0, 


13+0. 


01 


0.037+0.006 


0.93+0, 


04 


0.64+0 


,04 


1.4+0.1 


3.1+0.2 


0.14+0.01 


IRAS 08242-5050 


7. 


.68+0 


06 


0, 


.17+0. 


16 


0.098+0.063 


2.9+1, 


8 


1.9+0 


,8 


3.9+1.8 


5.0+2.4 


0.40+0.68 


IRAS 12553-7651 


7. 


.66+0. 


02 


0. 


18+0. 


04 


0.026+0.005 


0.6+0, 


1 


0.5+0, 


,1 


1.1+0.1 


3.8+0.9 


0.13+0.03 


IRAS 15398-3359 


7. 


.68+0. 


,01 


0. 


13+0. 


01 


0.22+0.03 


5.1+0, 


6 


4.2+0, 


,8 


8.8+1.7 


6.0+2.0 


0.20+0.06 


IRS 42 


7 


.72+0. 


,01 


0, 


12+0. 


03 


0.041+0.012 


0.82+0, 


11 


0.62+0 


,22 


1.4+0.4 


7.7+2.4 


0.21+0.06 


IRS 43 


7. 


.70+0 


,01 


0, 


21+0. 


04 


0.050+0.014 


1.7+0, 


6 


0.94+0, 


33 


2.1+0.6 


6.6+2.3 


0.70+0.36 


IRS 44 


7. 


,71+0. 


,01 


0, 


19+0. 


01 


0.053+0.015 


1.7+0, 


4 


0.9+0, 


,2 


2.1+0.4 


6.1+1.4 


0.63+0.23 


IRS 63 


7 


.70+0 


14 


0, 


.14+0. 


13 


0.042+0.026 


1.0+0, 


4 


0.72+0, 


,31 


1.6+0.7 


7.9+3.8 


0.42+0.59 


GSS 30 


7. 


.72+0. 


16 


0. 


14+0. 


13 


0.045+0.027 


1.1+0, 


.8 


0.7+0, 


,4 


1.6+0.9 


11+6.2 


0.60+ 0.9 


VSSG17 


7. 


.68+0. 


01 


0. 


.26+0. 


04 


0.042+0.006 


1.8+0, 


,2 


0.8+0, 


,1 


1.8+0.2 


11+2.2 


1.5+0.35 


RN091 


7 


.69+0 


,03 


0, 


27+0. 


13 


0.038+0.009 


1.6+0, 


7 


0.68+0 


,17 


1.6+0.5 


4.8+1.6 


0.76+0.44 


B59 YS05 


7. 


.70+0. 


.01 


0, 


19+0. 


06 


0.096+0.050 


3.2+3, 


4 


1.8+0, 


,9 


3.8+2.0 


2.7+1.6 


0.31+0.37 


2MASSJ17112317 


7. 


.69+0. 


01 


0, 


17+0, 


06 


0.13+0.052 


3.9+2, 





2.5+1, 





5.4+2.2 


2.8+1.2 


0.20+0.23 


SVS4-5 


7 


.67+0 


01 


0, 


16+0, 


02 


0.083+0.010 


2.2+0, 


2 


1.6+0 


,3 


3.5+0.6 


6.1+1.7 


0.29+0.08 


EC92 


7. 


,71+0. 


01 


0, 


07+0, 


10 


0.066+0.032 


0.67+0, 


21 


0.88+0 


.47 


2.0+1.2 


13+7.7 


0.33+0.23 


OA IRS5 


7. 


.68+0. 


01 


0. 


19+0, 


02 


0.055+0.011 


1.7+0, 


3 


1.1+0, 


,1 


2.3+0.3 


6.2+1.0 


0.50+0.13 


CrA IRS7b 


7 


,69+0 


.02 


0, 


16+0, 


05 


0.080+0.019 


2.2+0, 


5 


1.5+0 


5 


3.1+1.1 


3.0+1.2 


0.19+0.08 


OA IRAS32 


7. 


,72+0 


.04 


0, 


.09+0, 


13 


0.093+0.018 


1.4+0, 


4 


1.6+0 


,2 


3.5+0.4 


6.6+2.5 


0.12+0.05 


L1014 IRS 


7 


68+0. 


01 





19+0, 


06 


0.11+0.022 


3.7+1, 


.2 


2.4+0, 


,7 


5.1+1.5 


7.1+2.3 


0.53+0.25 


Average 




7.69 






0.16 




0.10 


2.7 




1.5 




4.1 


5.8 


0.34 


Standard deviation 




0.03 






0.05 




0.09 


2.0 




1.4 




3.5 


2.7 


0.35 
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3. Results 

3.1. CH 4 column densities 

Figure [T] shows the flux calibrated spectra of the 25 low mass YSOs containing the CH 4 
7.7 /im absorption feature. The spectra were converted to optical depth scale using a smooth 
spline continuum fitted to the 7-7.4 /im and 8.0-8.3 /im regions as shown in Fig. [IJ A local 
continuum is necessary since the CH4 feature lies on the edge of other broader features such 
as the 9.7 /im feature. Fitting the continuum is complicated by the presence of a weak 
HI recombination emission line at 7.5 /im (Pfund a) toward several of the sources e.g. the 
CrA sources, IRAS 03235+3004, L1489 IRS, RNO 91 and SVS 4-5. Of these ground based 
spectra around 4/xm exist for CrA IRS5, CrA IRS7b, IRAS 03235+3004, L1489 IRS and 
SVS 4-5. In all of these spectra the 4.05 /im HI Brackett a is clearly visible corroborating 
the identification of the emission feature with hydrogen. 

To estimate the uncertainties introduced by the choice of continuum, the spectra were 
also converted to optical depth scale by adopting a local, straight line continuum between 
7.4 and 7.9 /im. The differences in optical depth using the two different continua are less 
than 20% for most sources, but between 30 and 60% for a few cases: CrA IRS 7B, IRAS 
03235+3004, IRAS 12553-7651 and RNO 91. The resulting optical depth spectra using the 
smooth spline continuum subtraction are shown in Fig. [2l 

The peak position, full width half maximum (FWHM) and peak optical depth of the 
7.7 /im absorption feature along each line of sight were calculated from the best Gaussian 
fit between 7.5 and 7.9 /im (Table [2]). The reported error estimates include uncertainties in 
the fit and choice of continuum. The typical peak position is 7.69 /im and the peak widths 
range from 0.07 to 0.27 /im. The features hence range from unresolved to barely resolved 
as the resolution at the center peak position is 0.12 /on. The large widths of the feature in 
the observed spectra, except for Bl-b, suggest that CH 4 is generally in an H 2 or CH 3 OH 
dominate d mixture, where la boratory data show that the CH 4 feature has a width of up to 



0.15 /im (IBoogert et alJll997[ ). 



Figure [3] compares the source with the deepest absorption feature, Bl-c, with different 
CH 4 containing laboratory ice spectra, convolved to the resolution of IRS-SL1. In the three 
laboratory ice spectra used here - a H2O dominated ice, an ice mixture that contains equal 
parts of H2O, CH3OH and CO2 and pure CH 4 , all at 10 K - the CH 4 absorption profile has 
both different widths and peak positions. Toward most sources, as for Bl-c, H 2 dominated 
ice spectra provide the best fit. Hence, comparisons with a H 2 0:CH 4 3:1 ice spectrum were 
used to determine the amount of CH 4 present in the observed spectra and also how much of 
the 7.7 /im feature can be accounted for by solid CH 4 (Fig. [2]). Table [2] shows that CH 4 can 
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account for 45-100% of the absorption. 

The CH 4 column densities were calculated from the integrated optical depths of the 
laboratory spectra, scaled to the peak optical depths of the observations, and the band 
strength for the ben ding mode of solid CH4 in a H2O rich ice, 4.7 x 10 -18 cm molecule -1 



(IBoogert et al.lll997l ). The uncertainty in CH4 column densities stems from both the baseline 
subtraction and the uncertainty in ice mixture composition even after a H2O rich ice has been 
assumed. To obtain abundances with respect to solid H 2 0, the CH 4 column densities were 
divided by the solid H2O column densities from Paper I. FigureH] shows the CH4 abundances 
with respect to H2O as a function of H2O column density. The CH4 abundances in the entire 
sample vary between 2 and 13% and it is seen in the plot that the sample can be split into 
two parts: sources with H 2 column densities around 2xl0 18 cm -2 and sources with H 2 
column densities of 3-40 xlO 18 cm -2 . In the former group the CH4 abundances with respect 
to H 2 vary between 6 and 13% and in the latter group all CH4 abundances fall between 2 
and 8%. Due to the low total column densities in the high abundance group, the uncertainties 
there tend to be larger. Below 3xl0 18 cm -2 there also seem to be some negative correlation 
between column density and CH4 abundance. FigureH] also shows the CH3OH abundances 
and upper limits towards the sources in this sample (Paper I) . CH 3 OH abundances span a 
larger interval than CH 4 and show none of the column density dependences visible for CH 4 . 
The large variation for CH3OH abundances is similar to what was found previously for a 



different sample by IPontoppidan et al.l (120031 ) . 



The excess absorption (0-55%) of the observed 7.7 fim features in many of the astro- 
nomical objects is due to broader profiles than expected even for CH4 in a H 2 rich ice. It 
is possible that an additional molecule is contribu ting to the op t ical d epth of the 7.7 jttm 



feature, e.g. solid S0 2 . Solid S0 2 was suggested by lBoogert et al.l (Il997f l to explain the blue 
wing of the 7.7 /im feature in W33A. In contrast, another high mass source NGC7538 IRS9 
displays no such wing. A comparison between these two sources and low mass sources from 
this sample shows that the same variation is present here (Fig. [5]) - approximately 25% of 
the sources in this study have a clear blue wing, perhaps attributable to solid S0 2 . The max- 
imum amounts of S0 2 present in the observed ices is estimated by assuming that all excess 
absorption is due to SQ 2 and us ing its measured band strength of 3.4 x 10 -17 cm molecule -1 



dSandford fc Allamandolal Il993h Table [2]). The solid S0 2 abundances then vary between 



0.1 and 1.5% with respect to H 2 0. 

In cases where the excess is similar on both sides of the CH 4 absorption, another source 
of the large widths of the features may be the choice of continuum; the spline continuum 
was difficult to fit because the CH4 feature is generally shallow and overlapping with other, 
larger, features. It is possible that the wings toward several sources are not intrinsic, but 



a product of this fit. This is highlighted by a comparison between Fig. [I] and Fig [21 which 
shows that the sources with clear continua also have thinner absorption features than the 
average. 



3.2. Upper limits of solid CH 4 

Upper limits for solid CH 4 were determined for the 27 sources without CH 4 detections, in 
the sample of 52 low mass ice sources originally probed for CH4, by estimating the maximum 
amount of a H 2 0:CH 4 3:1 ice that could hide under the noise. The average 3cr upper limit 
is 15% solid CH 4 with respect to solid H 2 0, which is in the upper range of CH 4 abundances 
in the sources with solid CH 4 detections. All abundance upper limits below 30% are shown 
in Figure 4. Only one CH 4 upper limit, of ~3%, falls below the average abundance of 4.8%. 
The lack of detection in these 27 sources is hence probably due to the spectral quality and 
low fluxes of the objects rather than a lower amount of solid CH 4 . 



3.3. Molecular correlations 

The objects in this sample may vary significantly in temperature structure as well as 
other environmental factors. Lack of correlations between molecular abundances may hence 
be due to either different formation pathways or differences in volatility. A lack of correlation 
between CH 4 and molecules of different volatility does hence not exc lude that they forme d in 



a similar manner. In pure form CH 4 has a similar volatility to CO (jCollings et al.l 120041 ) . If 
CH 4 is mixed with H2O significant amounts can be trapped inside of the H2O ice, however, 
and then CH 4 has an effective volatility closer to that of H2O and CO2. 

The column densities of solid H 2 (26 detections) and CH 3 OH (10 detections) have been 
derived in Paper I a nd so lid CO (13 ground based observations) and CO2 (25 detections) in 



Pontoppidan et al.l (120031 ) and Paper II for many of the CH 4 sources. The column densities 
and abundances for these molecules are plotted versus CH 4 in Figs. [6] and [71 Where the 
plots reveal a linear relationship between molecules the best linear fit is also drawn. The 
Pearson product-moment correlation coefficient, R, was calculated to measure the strength 
of the correlation - R 2 directly gives the fraction of variance of the second molecule that is 
due to changes in CH 4 , assuming a linear relationship between the two molecules. 

Figure M shows the correlations between the column densities of solid CH 4 and the 
four other ice components. Some correlation between column densities of ice species and 
with total column density is always expected, but Fig. [6] shows that the strength of this 
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correlation is variable between CH4 and the different molecules. CH4 correlates strongly 
(_R 2 =0.91) with CO2, which is believed to form on grain surfaces. The correlation with H 2 0, 
another species formed on surfaces, is weaker (i? 2 =0.64), but this correlation coefficient is 
significantly increased if one outlier, IRAS 03245+3002, is removed. R 2 is then 0.82 for 
the CH4-H2O correlation. IRAS 03245+3002 could have been identified previous to the 
correlation studies as an outlier since it is the only source with no CO2 ice detection putting 
an upper limit on its C0 2 /H 2 of 15%, compared to 30+9% for the entire sample. Even 
without removing the outlier the correlations are significant at the 99% level with 23 and 24 
degrees of freedom (DF) respectively. 

CH 4 and CH 3 OH and CO are barely signficantly correlated (i? 2 =0.48 with 8 DF and 
0.49 with 11 DF respectively). The small number of CO abundances, which are only available 
for 11 of the targets in this study, complicates interpreting the low correlation plots between 
CO and CH4. If only the sources for which CO measurements exist are used to calculate the 
CH4-CO2 and CH 4 -H 2 correlations these are reduced to i? 2 =0.46 and 0.74 respectively. 

The stronger correlation between CH 4 and C0 2 compared to CH 4 and H 2 is curious, 
but may be an artefact of the fact that the C0 2 column densities are better known for this 
sample than the H 2 column densities. For most sources in this sample the H 2 column 
densities are determined by estimating the depth of the H 2 libration feature, which is 
difficult due to the overlap with the silicate feature at 9.7 /im. The uncertainties in these 
H 2 column densities were estimated to ~30% in Paper I after comparison between derived 
abundances from the 3 /zm mode and the libration mode for the sources with 3 /zm spectra. 
Furthermore the different amounts of C0 2 toward the different s ources may in t roduc e further 



error in the derived H 2 column densities as recently shown in lOberg et al.l (120071 ). 



Figures [7J show the correlations between solid CH 4 and the other three ice species, with 
the abundances normalized to the solid H 2 column density - the main ice constituent (Fig. 
[7]). When normalizing to H 2 there is no significant correlation at a 95% level between CH 4 
and CH3OH and CO (R 2 below 0.1) and at best a weak correlation between CH4 and C0 2 
(R 2 = 0.27). There is also no correlation between CH4 abundances and any of the three 



CO components as defined by iPontoppidan et al.l (120031 ). with the different components 



corresponding to pure CO and CO in a H 2 and C0 2 rich ice respectively. 



3.4. Spatial trends 



The average CH4 abundance relative to H 2 and its standard deviation is shown for 
all clouds with more than one detection and for the entire sample in Table [3J The averages 
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for the individual clouds range from 2.8% to 9%, with Serpens and Ophiucuhs at 9 and 
8% respectively. Whether the clouds are significantly different in their CH 4 abundances was 
evaluated using Analysis of Variance with the StatisticslOl resampling software. Resampling 
is more robust than traditional statistical tests since there is no need to assume an underlying 
distribution, which is especially useful when the sample size is small. The test procedure 
starts with calculating the sum of the absolute deviations of the cloud averages from the 
sample average. The CH 4 abundances in the five clouds are then randomly resampled into 
five new groups with the same group size distribution as before and the absolute deviations 
of the group averages from the sample average is re-calculated. This is repeated 1000 times 
and the number of times the resampled sum of deviations exceeds the sum of deviations of 
the cloud averages is counted. The nul-hypothesis that the difference between the clouds is 
due to chance could then be rejected with 95% confidence. If Ophiuchus is removed from 
the sample there are no longer any significant differences between the clouds. 
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Table 3. The average CH4 abundance relative to H2O, and its standard deviation for five 
of the investigated clouds and for the entire sample. 



Cloud 


Number of sources 


Average 


St. dev. 


Perseus 


6 


0.040 


0.014 


Ophiuchus 


6 


0.083 


0.019 


B59 


2 


0.027 




CrA 


3 


0.052 


0.019 


Serpens 


2 


0.093 




All 


19 


0.058 


0.027 
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Discussion 



4.1. Low vs. high mass YSOs 



For comparison the solid CH 4 abundances with respect to H 2 for several high mass 



YSOs (Table [T]) were re-derived and plotted relative to H 2 abundances from iGibb et al. 



(2004) in F ig. [H These source s have been previously investigated for solid CH 4 by lBoogert et al 
(119961 ) and lGibb et al.l (120041) . The so l id CH 4 abundances derived here are 30 to 60% higher 
than those published by Gibb et al. _( 200-1) ma inly due to a difference in band strength 
(7.3xl0~ 18 molecules cm -1 by lGibb et al.l (12004 ) and 4.7xl0 -18 molecules cm -1 here). The 
plot shows that the high mass ISO CH 4 sources fit in seamlessly with the low mass CH 4 
sources with H 2 column densities above 2xl0 18 cm -2 . Figure [5] also shows that there is 
a similar variation in the 7.7 fim feature profile between low and high mass objects. There 
is hence no reason to expect large systematic differences in CH 4 abundances or formation 
pathways toward low and high mass YSOs. 



4.2. Formation scenarios 



The similar or higher abundance of CH 4 with respect to H 2 toward the low mass sources 
in this sample compared to what has been found previously toward high mass young stellar 
objects, suggests that the CH 4 formation rate is not dependent on stellar UV irradiation. Of 
the three formation scenarios for CH 4 suggested in the introduction, this study hence does not 
support a formation pathway connected to stellar UV processing of CH3OH. Furthermore, in 
the sources with 2-40 xlO 18 cm~ 2 H 2 the CH 4 abundance is nearly constant. In comparison 
the CH3OH abundances vary by a factor of 10 and it seems unlikely that this would be the 
outcome if the main formation pathway of CH 4 is connected to CH 3 OH through e.g. cosmic- 
ray induced UV processing of CH 3 OH, which is also present under quiescent conditions. 

Of the two remaining scenarios, formation in the gas phase with subsequent freeze- 
out and hydrogenation of C on grain surfaces, this study lends most support to the latter 
mechanism. First, under quiescent conditions, gas phase models predic t steady state total 
CH 4 / H 2 abundances of only around 10~ 7 — 5 x 10™ 14 (IWoodall et al.l 120071 ; iBergin et al. 
Il995[ ). compared with our inferred CH 4 /H 2 abundances of ~ 2 — 13 x 10 6 (assuming a 
standard H oO/H? ratio of 10 ~ 4 ). In early times when C/CO>l, the CH 4 /H 2 can reach 
above 10~ 6 (IMillar et al.lll99lT) . Ices cannot form at extinctions lower than ~ 2 Ay, however 
( Cuppen fc Herbstl 2007 ). At extinctions higher than 2 Ay, C/CO<l and hence freeze-out of 



the high CH 4 abundances at early times cannot be used to explain the high i ce abundances. 



In addition, pure CH 4 has a similar volatility to CO, within a few degrees (jCollings et al. 
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2004j ). and if both are formed in the gas phase and subsequently frozen out, the two molecules 
should correlate, both in absolute column densities and in abundances relative to the much 
less volatile H 2 0. This study clearly shows that this is not the case. 

This weak correlation between CO and CH4 can be contrasted with the stronger corre- 
lations between CH4 and CO2 and H2O column densities. Furthermore the lack of significant 
correlation with any species once the correlation with H 2 has been effectively divided out, 
by normalizing with the H 2 column densities, shows that CH 4 is not significantly better 
related to any other molecule than to H 2 0. Together this suggests a formation scenario of 
CH4 more related to H 2 and C0 2 formation, which both form on grains, rather than to CO. 
In theory the correlations could be single-handedly due to the fact that CH4 has a higher 
desorption temperature than in its pure form due to mixing with H 2 0. This, however, is only 
expected to occur if CH 4 is formed together with H 2 and hence CH 4 formation on grain 
surfaces is still the most plausible outcome of this study. 

The hydrogenation scenario is also supported by the broad profiles of the features which 
have only been found in laboratory spectra when CH 4 is in a mixture with hydrogen bonding 
molecules, which are believed to form on grain surfaces. These conclusions are in agreement 
with studies of the solid CH 4 profiles t oward high mass YS Os, which show that solid CH 4 



is mixed with H 2 or possibly CH3OH (IBoogert et al.lll997l) a nd also with the observed low 



gas/solid ratio of CH 4 compared to CO (IBoogert et al.ll2004al ) 



4.3. Differences between clouds 

Figure @] and Table [3] shows that the CH 4 abundances differ from cloud to cloud. In Fig 
HI 4 of the 5 sources with high CH 4 abundances and H 2 column densities below 2xl0 18 
cm -2 are Ophiuchus sources (the fifth is in Serpens). The two Ophiuchus sources that have 
higher H 2 ice abundances also have more 'normal' CH 4 abundances. This indicates that 
it is the low total column density towards four of the Ophiuchus sources rather than initial 
conditions in Ophiuchus as a whole that is responsible for the extreme CH 4 abundances 
found there. This may be due to that towards sources with low total column densities the 
C/CO ratio may be larger allowing for more CH 4 to form. The C/CO ratio is dependent on 
extinction for low extinctions, but not for high ones. This is consistent with our observation 
that CH 4 abundances are only dependent on total column densities for low column densities. 
Sakai et al. (ApJ in press) have instead suggested that different CH 4 abundances could be 
due to different collapse times. In clouds that collapse very fast, chemical equilibrium may 
not be reached favoring a C-based rather than CO-based chemistry. This question can only 
be settled if gas phase data for these sources become available. 
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4.4. Comparison with models 



S olid CH4 abundance s have been modelled previously for a variety of conditions (IHasegawa et al. 



19921 ; lAikawa et al.ll2005l ). The predicted CH 4 ice abundances with respect to H 2 usually 
vary between psI-10%, with the main formation path being sequential hydrogenation of C 
atoms on grains at a time when a large fraction of the gas phase C has already been con- 
verted to CO, i.e. C/C0<1. The observed abundances in our sample fall mostly within this 
range and are hence in general agreement with the models. 



Aikawa et al.l (120051 ) model molecular abundances on grain surfaces and in the gas phase 



during the collapse of a spherical cloud. Regardless of initial conditions and the central den- 
sity when the collapse is stopped, the CH4/H2O and CH4/CO2 ice ratios are fairly constant 
as a function of radius in the cloud. The value of CH 4 /H 2 ratio varies with initial condi- 
ti ons between a few and 1 0%. In contrast the CH4/CO ratio is radius dependent. The model 
of IHasegawa et al.l (119921 ) is more aimed at modelling interstellar clouds before collapse, but 
as they show the time evolution of the chemistry it may still be useful to compare with our 
observations. As in the collapse model, CH 4 , C0 2 and H 2 trace each other fairly well, while 
CO and CH4 are only correlated during some time intervals, regardless of initial conditions. 
Both these models are consistent with Fig. [6] where solid CH 4 is better correlated with H 2 
and C0 2 than with CO. 



5. Conclusions 

We present Spitzer-IRS spectra of the solid CH 4 feature at 7.7 /zm toward a large sample 
of low mass young stellar objects. Our conclusions are as follows: 

• 25 out of 52 low mass young stellar objects show a solid CH 4 feature at 7.7 /xm. 

• The solid CH 4 abundance with respect to H 2 is centered at 5.8% with a standard 
deviation of 2.7% in the sources with CH 4 detections. In the sources without detections 
the average upper limit is 15%, which is not significant compared with the detections. 

• The sources (two Ophiuchus and one Serpens) with more than 10% CH 4 all have H 2 
column densities below 2xl0 18 cm~ 2 . Due to the low total column densities, two 
of these three sources have uncertainities larger than 50%. Above 2xl0 18 cm -2 the 
sources (20 out of 25) have a fairly constant CH 4 abundance of 4.7±1.6%. 

• The 7.7/im feature profiles are signficantly broader for all but one object than expected 
for pure solid CH 4 and toward most sources also broader than expected for CH 4 in H 2 
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dominated ices. Approximately 30% of the features have a blue wing, seen previously 
toward high mass YSOs and there attributed to solid S0 2 

• The column densities of solid CH 4 and H 2 and C0 2 are clearly correlated, while CH 4 
and CO and CH 3 OH are only weakly correlated. 

• There is also no correlation between the CH 4 and CO abundances when both have 
been normalized to the H 2 abundance. 

• The Ophiuchus cloud has significantly higher CH 4 abundances compared to the rest 
of the sample, probably due to the low total column densities towards several of the 
sources. There are no significant differences between the remaining clouds. 

• The abundance variation is smaller for CH 4 compared to solid CH 3 OH; CH 4 seems to 
belong to the class of molecules, also including H 2 and C0 2 that appear 'quiescent', 
i.e. their abundances are more or less constant, in contrast to highly variable ices like 
CH3OH and OCN~. If the Ophiuchus sources are included CH 4 is somewhere between 
the two classes. 

• Sample statistics and comparison with model predictions support CH 4 formation through 
hydrogenation of C on grain surfaces. 
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Fig. 1. — The Spitzer-IRS spectra (solid line) of the CH 4 sources between 7.0 and 8.2 /xm, 
plotted together with the chosen spline continua (dotted line). The tick bar in the upper left 
corner indicates the noise level derived within the c2d pipeline. 
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Fig. 2.— Optical depth spectra of the CH 4 sources plotted together with 
laboratory spectra (dashed) of a CH 4 :H 2 mixture 1:3 (Leiden databases at 
www.strw.leidenuniv.nl/~lab/databases) convolved to at the same resolution and sampled 
in the same way as the observational data. 
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Fig. 3. — Optical depth spectra of Bl-c (solid line) plotted together with laboratory spectra 
of a CH 4 :H 2 1:3 mixture (dashed), a H 2 0:CH 3 OH:C0 2 :CH4 0.6:0.7:1:0.1 mixture (dotted) 
and pure CH 4 ice (dashed-dotted). The derived column densities using the three different 
ice compositions are 1.7, 1.6 and 1.7xl0 18 cm" 2 respectively. 
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Fig. 4. — The CH 4 and CH 3 0H abundances relative to H 2 are plotted versus the column 
densities of H 2 for Spitzer IRS (black diamonds) and ISO (grey stars) CH 4 sources. Two 
groupings are visible in the upper plot: sources with H 2 column densities around 2xl0 18 
cm -2 with high CH 4 abundances and sources with H 2 column densities of 3-40 xlO 18 cm -2 
and a nearly constant CH4 abundance around 4-5%. CH3OH abundances show no similar 
groupings and span a larger interval than CH 4 . CH 4 and CH3OH abundance upper limits 
below 30% are plotted in grey. 
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Fig. 5. — 7.7 fim profile comparison between high mass sources W33A and NGC7538 IRS9 
(with their original SWS resolution and convolved to the resolution of IRS-SL1) and four low 
mass sources. The two Gaussians plotted together with W33A fits the observed spectrum, 
but only the thin component centered at 7.7 /im is explained by CH 4 . Several of the low 
mass sources have similar profiles to W33A with an additional component, attributed to 
S0 2 , that makes up the blue wing of the 7.7 /im feature. A few low mass sources also show 
thin profiles similar to that of NGC7538 IRS9. 
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Fig. 6. — The column densities of four ice species, CH 3 OH, H 2 0, CO and C0 2 are plotted 
versus the column densities of solid CH 4 . Solid CH 4 is strongly correlated with solid C0 2 
and H 2 0. CH 4 is only weakly correlated with solid CO and CH 3 OH, but this may be partly 
due to the fact that CO has only been observed and CH 3 OH only detected towards much 
fewer targets. The outlier in the H 2 plot is IRAS 03245+3002. 
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Fig. 7. — The abundances of solid CH 3 OH, CO and C0 2 relative to H 2 are plotted versus 
those of solid CH 4 . There is no correlation between the relative amounts of solid CH 4 and 
CH3OH or CO. In contrast the CH 4 and solid C0 2 abundances are weakly correlated. 



